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INTRODUCTION 


equations of the form 


The dynamic response of v/1nd turbines to turbulent 
wind fluctuations has generally been modeled using 
only changes in the streamwise wind velocity as the 
disturbance which causes the varying aerodynamic 
forces and moments. Often these wind fluctuations 
are thou^it of, and modeled, as discrete wind 
gusts of specified shape and duration, which occur 
at some average rate. In addition, these wind 
fluctuations usually have been assumed to act uni- 
formly over the entire rotor disk. In terms of 
the designer's needs, these models are used to 
generate design loads and expected control system 
actions. 


It was the objective of the work reported here, 
and in the companion paper [ 1 ], to take a broader 
look at wind turbine dynamic response to turbu- 
lence, and attempt to ascertain the features of 
turbuience that wind turbines are most sensitive 
to. A statistical description of the wind input 
including all three wind components and allowing 

linear wind gradients across the rotor disk, was 

used together with quasi-static aerodynamic theory 
and an elementary structural model involving only 
a few degrees of freedom. The idea was to keep 

the turbine model simple and show the benefits of 

3 ::.-:.:^..-.:. this type of statistical wind representation be- 
fore attempting to use a more complex turbine 
model. As^ far as possible, the analysis was kept 
in the simplest form, while still preserving key 
physical responses. 

From the onset of this work, it was felt that the 
results should be validated by comparison with test 
measurements . Due to the three-bladed rigid rotor 
used on the turbine and the limited degrees of 
freedom, comparison with data from one of the small 
systems under test at Rocky Flats would provide the 
most realistic comparison. At this time, the ex- 
perimental comparison is incomplete. 

The Turbine Model 


The" wind turbine model is shown schematically in 
Figure 1. Both the rotor and the nacelle are 
assumed to be rigid bodies which move in unison, 
except .for_Lhe spinning rotor. Due to tower 
flexibnity, the nacelle and rotor are free to 
translate in a plane parallel to the ground and 
rotate about the top of the tower in pitch and 
yaw. The yaw angle of the rotor axis 1s defined 
by the angle, <|>, and the pitch angle by x- The 
lateral translation, U, is in the x direction, 
while the V translation is in the y direction 
along the rotor axis. The rotor spin velocity 
is given by Q + f , where q is the mean rotation 
rate and f is some small fluctuation. For the 
casr 3 of a turbine with a three-bladed rigid rotor, 
the basic principles of Newtonian mechanics and 
linear, quasi-steady aerodynamics give motion 


K.jXj * C,.t. . K,jXj 


F. .V. 
J 


(i: 


where Mij, Cjj, and F-jj are the turbine system 
inertia, damping, stiffness and wind input coeffi- 


cients, The five displacement coordinates already 
described are X j , while the wind inputs are Vj . 


The Tower 


The wind turbine tower was modeled as a single 
finite element within which the tower displace- 
ments were expressed in terms of interpolating 
polynomials and the displacements at the top of 
the tower. Then the tower deformation, v(z,t), 
about one bending axis was written in the form 

v(z,t) = Py(z)V(t) + P^(z)x(t) (2) 

where Py and Px are the interpolating functions 
which approximate the displacements within the 
tower. These are conveniently represented as 
cubic polynomials satisfying the necessary bounda- 
ry conditions of a cantilever tower. Using this 
expression for the tower bending displacement, the 
stiffness and mass coefficients may be obtained by 
one of the numerous energy methods. In terms of 
the interpolating functions Py and the gener- 
alized stiffness and mass coefficients for the 
tower may be expressed as 
L 


k.. = f EI(z)PV(z)PV(z)dz 

ij 0 ^ ^ 

(3) 

L 


m.. = / m(z)P.(z)P.(z)dz 
0 ^ 3 

(4) 


where EI(z) and m(z) are the stiffness and mass 
per unit length as a function of height. For ad- 
ditional detail concerning this technique, the 
interested reader should see Clough and Penzien 
[2]. Although the tower properties are the same 
in both bending directions, only one degree of 
freedom was desired for the x direction and there- 
fore rotation of the nacelle about the rotor spin 
axis was neglected. The method of static 
condensation was used to eliminate the unwanted 
degree of freedom and obtain the desired x direc- 
tion stiffness and inertia coefficients as 


kyu = '^vv ^ 


(5) 


'"uu " ""vv ■ 


In addition, the nacelle and rotor inertias add 
directly to the tower inertia coefficients, mjj, 
to give the turbine system inertias. A detailed 
tabulation of the various terms in the inertia 
matrix is provided in the Appendix. There is also 
a gyroscopic coupling between the turbine 
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pitch and yaw motions. This coupling coefficient 
appears in the damping matrix and is given by 


<(>X 


= 


-C 


x<l> 


(7) 


where I is the total effective inertia of all the 
spinning mass connected to the turbine rotor. 

Using this simple model the structural stiffness 
and Inertia coefficients for a particular wind 
turbine are numerically calculated using a TI-59 
calculator code [3]. 


Aerodynamic Forces 


The geometry of the three-bladed rigid rotor is 
illustrated in Figure 2. The blades are coned at 
an angle Bq, and are assumed to be twisted and 
tapered. The angle 6 defines the pitch setting 
as the angle from the plane of rotation to the 
zero-1 ift-line of the airfoil at the blade tip. 


For this analysis, quasi-static aerodynamics will 
be used to compute the forces acting on the blades 
due to the turbulent wind and structural motion. 
The wind input including turbulence is assumed to 
be made up of a steady mean wind, V^, plus fluctu- 
ating components, V-j(t), which at any instant are 
constant over the rotor disk and turbulence gradi- 
ents Vt i(t), which vary linearly across the disk. 
Both Vi(t) and Vi^j{t) may be thought as disk 
averaged time dependent quantities. This allows 
the wind velocity to be written in a linear expan- 
sion as 






^x,y ^x,z' 
jjy^y }y,z 
^z,y Vz,z, 



(8) 


where r is the radial position in the rotor disk, 
and Qt is the azimuthal location. Motivation for 
this particular form for the turbulence and dis- 
cussion of its accuracy is presented in the com- 
panion paper, reference [1]. In this equation 
the mean wind direction coincides with the y axis 
of Figure 2. In addition, the spatial change in 
the y direction turbulence component due to the 
coning has been dropped In the above expression, 
which eliminates the effect of the turbulence 
gradients Vi ,y. 


+ sinl^t[-x(V^-v.) + 6V^ 

- ax]}/Rn 

= (V^-v.)/RQ 4- {6Vy+ cosrit[-BQ6V^ + rx] 

+ sinQt[3Q(U-5V^) 

- r4)]}/Rf^ (10) 

where the fluctuating part of the wind turbulence 
has been written as 6V-j to shorten the expressions. 
These expressions have also been linearized assum- 
ing small displacements, and wind fluctuations; 
however, in some places the product of the static 
coning angle and the wind fluctuation were retained 
because of interest in their effect. 

Referring to Figure 3, the aerodynamic forces 
parallel and perpendicular to the rotor plane may 
be written a$ 


dA^ = 1 pa- c d?{-nvj - ev^v^ + (l-ri/2)v2} 

( 11 ) 

dA^ = 1 pa- c dC{-9V^ + 

where the lift and drag for a blade element have 
been calculated using the static formulas with the 
instantaneous velocities. In the above expressions 
a’ is the slope of the sectional lift curve 
dCt/da, c is the local airfoil chord, n is the 
ratio Co^/a', and p is the air density. 

Using the wind input of Eq. (8), together with the 
velocity expressions of Eq. (10) and substituting 
into the aerodynamic force relationships gives 

dA^ = j pa 'c^(Rfi)^Rdx{A’-B‘ [f^+f^cosQt+f^sinf^t 

+fc2C0s2^^t+fs2Sln2fit] 

+C‘ [g^+g^cos^zt+g^s1n^2t 

dA^ = pa ’c^(R^^)^Rdx{D*+E' [f^+f^cosQt+f^sinnt 
+fc2COs2Qt+fs2Sin2r6t] 

+F'[g +g cosQt+q sinQt 

‘-^0 ^C 'S 


Using the above wind representation and funda- 
mental kinematic relationships provides the rela- 
tive velocity as observed from the turbine blade. 
The relative velocity is made up of contributions 
from the wind, the moving structure and the in- 
duced velocity caused by aerodynamic action. In 
equation form this is 



In terms of the displacement coordinates, the rela- 
tive velocity components parallel and perpendicular 
to the rotor disk rotational plane are given by 

y^/RQ = rJl/RQ +{r'? + cosf^t[-4s(V^-v^ ) - 6V^ 

+ 0 - a^.] 


+g^2^os2nt+g^2S^^2^2t]} (12) 

where the primed quantities are the aerodynamic 
constants 

A'(x) = [(l-r]/2)A^-x(nx+eX)]c/c^ 

B'(x) = [2nx+0X]c/c^ C’(x) = [{2-n)X-ex]c/c^ 

D'(x) = [(1+n)Ax-e(x^+>X/2)]c/c,. 

E'(x) = [(l+n)x-20x]c/c^ 

F'(x) = [(1+n)x-xe]c/Cj. (13) 

with x = r/R = (h+c)/R. X = (V -v.)/Rn, and c. is 
the chord at the rotor tip. NljteHha.t both the 
pitch setting, e, and the blade chord, c, may be 
functions of x. In the above force equations the 
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subscr'^pted f and g variables are combinations of 
the wind Inputs and response variables and are de- 
fined as follows: 


f = -4^X+(-V + 0 - ai)/RQ 
c X 

^s2= 

',x- I«r.x - \.z> 

'zx" 2^^z,z ~ '^x,x^ 

9 o = (Vy - V - 


= {rx + rVy^^ - BgV^)/Rn 
9s = (6o(0 - \)-rl ^ rVy_^)/Rfl 


The second wake model was called the '’equilibrium 
wake," For this model, the axial fluctuations in 
wind velocity are assumed to occur so slowly that 
the induced velocity is the steady state value for 
the instantaneous wind speed. In this situation, 
the axial flow will be not only time varying but 
nonuniform, because of the inclusion of the fluctu- 
ating wind gradient terms in the. turbulence model 
of Eq, (8). These gradients could be thought of as 
slowly changing wind shears of arbitrary orienta- 
tion, since their effect on the wind turbine 
would be similar. To obtain an approximation for 
the induced velocities of this "equilibrium wake," 
the "semi-rigid" wake model discussed by Miller 
[5] was used. Miller shows that the effect of 
Including the induced velocity due to the nonuni- 
form flow is to reduce the lift by a factor re- 
ferrred to as the "lift deficiency" function. 

For this analysis, assuming small velocity changes, 
the lift deficiency function was approximately 

^d(x) " 1 +t^F'(x]/x(2A-A^J 


9 c 2= - ^2 


'zx 

^zx 


" ^cZ 

’ 2^^z,y. * '^x,z^ 

■ K.x * 


(14) 


To obtain the aerodynamic coefficients for the 
total forces acting on the rotor hub, the appropri- 
ate components of the blade element forces, Eq. 
(12), are summed over the three blades and inte- 
grated with respect to radius for a specified in- 
duced velocity distribution. This gives the net 
thrust, torque, horizontal and vertical forces, 
and the yaw and pitch moments, which are to be 
added to the structural terms resulting in the 
final equations of motion, Eq. (1). A detailed 
list cf these equations is provided in the 
Appendix. 


The In duc ed Velocity 

The”“aerodynarrti cs of wind turbines involve highly 
complex flow phenomena, which require rather 
sophisticated theories in order to obtain accurate 
predictions. However, some fairly simple theories 
making relatively crude assumptions can often give 
reasonable estimates and generally can give excel- 
lent Insight into the physical phenomena of inter- 
est, In this case two different wake models were 
used in an effort to gain Insight into the signi- 
ficance of changes In the induced velocity field 
on wind turbine response to turbulence. 

For the first wake model, the induced velocity was 
computed using blade element theory following the 
approach of Wilson [4], and performing a momentum 
balance neglecting wake rotation. This provided 
the induced velocity as a function of radius, 
under the assumption that the rotor axis was per- 
fectly aligned with the wind direction. After 
the induced velocity distribution was computed 
for a given mean operating condition, it was as- 
sumed to be constant and independent of turbu- 
lent wind fluctuations. This model was named 
the "frozen v^ake model." 


where X = V and t. = 3a’c./87rR. In addition, 
the aziJi(uthaY change in the inSuced velocity dis- 
tribution led to a change in the in-plane aerody- 
namic coefficients B'(x), where the change was 

AB'(x) = t^[B’(x)F‘(x)+C'(x)E‘(x)]/ 

{x{2X-X^)> (16) 

Because this change is small, it was tempting to 
neglect it, but all of the In-plane forces are 
small so it was retained. Finally, the wind 
fluctuations in the axial direction V , V and 
V are associated with a change In momeritum in 
tK^^streamwi se direction which, for the assump- 
tions of this wake model, change the equilibrium 
thrust. This added lift factor Is approximately 

Xg(x) = 1 + A/(2X-X^) (17) 

To incorporate these effects, the aerodynamic 
coefficients B'(x), C‘(x), E'(x) and F’(x) of 
Eq. (13) were modified in the following manjier to 
obtain the "equilibrium wake" coefficients: 

B;(x) =X^(x){B'(x)+AB'(x)} 

Cg(x) =^^(x)C'(x) 

E;(x) =Jfj(x)E-(x) 

Fg(x) =Jfj(x)E'(x) 

M 

F^y(x)= Jtj(x)3C^(x)C'(x) (18) 

where the two coefficients Cey and Egy are speci- 
ficdlly associated with the wind fluctuations Vy, 
Vy^z and Vy^x* 1^*^^ aerodynamic coefficients A'tx) 
and D‘(x) are related to the mean thrust and torque 
and are thus unaffected by wind fluctuations. 
Computationally, the influence of the wake model ^ 
can be observed by changing the primed aerodynamic 
coefficients in the blade element force relation- 
ships Eq. (12). 
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Although Doth of these wake models were useful in 
developing an understanding of the influence of 
the induced velocity distribution on machine 
response to turbulence, it is unclear whether 
either model accurately approximates the real 
distributions and future work is needed to 
evaluate the effects of unsteady wake aero- 
dynamics. 

State Space Equations 

The equations of motion Eq. (1) can be written in 
matrix form 

[M]{X}+[C]{^}+[K]{X} = {Q^}+[F]{u} (19) 

where 


{X} = “ displacement coordinate 

{Ql}^= (0,T,0,Q,0) = steady state 

{u}^ = (V ,V ,V ,V ,V ,Y ,e ,Y .e ) 

X y z y,x y,z 'xz’ r 'r xz' 

= wind inputs 

cos3Qt + sinSfit 


{x} = [A]{x} + [B]{w) 
{y} = [c]{x} 


where 


{x} = 


( 21 ) 


'{X} 

^ [A] = 

'[0] 

-I 

[I] 

-I 

[0] 

{X} 

*[M] [K] 

-[M][K] 

[M][F] 

,{u}J 


JO] 

[0] 

[A ] 
w-* - 


[B] = 


[ 0 ] 

[B ] 
w-* 


{y} My r= outputs 
Power! 

Jx} J 


[c] = 


[K] 


[ 0 ] 


(oo...nCg...O) 

[I] 


= response matrix 


Yp = sin3nt + Y^)( cos3fit 

Furthermore, it is possible to model each of the 
nine turbulence inputs using a set of stochastic 
differential equations of the form 

iu) = [A^]{ul + [B^Kwl (20) 

where the components of {w} are white noise of 
equal power spectral density, [By^] is the white 
noise input distribution matrix which Is diagonal. 
The [A^] matrix is diagonal except for two ele- 
ments 


[A ] 
w-* 


dyy 3Q 

-3Q agg 

a 99 


which arise as^a result of the sinSr^t and cos3Qt 
In the c and y wind inputs. A more detailed 
presentation or the wind input model is pre- 
sented in the companion paper [1]. 


Discarding the steady terms, it is convenient to 
transform the equations of motion given in Eq. 
(19) to the state space form, so that they are 
written as a set of first order equations similar 
to the turbulence inputs of Eq. (20), To further 
facilitate the computation of results, the state 
space form of Eq. (19) can be augmented with the 
turbulence inputs, Eq. (20), to form a single 
system of equations with white noise as the 
driving input. The five turbine displacements 
and their derivatives together with the nine 
turbulence Inputs will form the state vector for 
this augmented system. The governing equations 
can then be written 


With this formulation it is a relatively straight- 
forward numerical procedure, to determine the com- 
plex eigenvalues of the A matrix and then to com- 
pute the modal matrix, which is made up of the 
associated eigenvectors. The modal matrix can then 
be used to decouple the equations of motion, so 
that transfer functions between any of the nine 
white noise inputs and any output, y-j , may be 
easily computed. These transfer functions account 
for differences in the energy level for the turbu- 
lence inputs, {u"^, so that a comparison of the 
transfer function magnitudes provides a direct 
estimate of relative importance. The final re- 
sult uses the central equation from random 
vibration theory, which states that the spectral 
density for any of the outputs {y} will be given 
by 

{Sy(o))} = [|Hy^(..)|2] {S^(O))] 

for uncorrelated inputs, where {Sy(w)} is the p 
spectral density of the outputs {y}, [[Hyw(aj)r] 
is the matrix of the transfer functions magnitude 
squared and {S^(w)} is the spectral density of the 
white noise driving inputs, which are all equal, 

Resul ts 


To determine the influence of the turbulence in- 
puts modeled in this work two wind systems of 
vastly different size were examined. The smaller 
machine, called the Mod-M, was an 8 kW turbine, 
with a three bladed rotor located downwind of the 
tower, and designed for free-yaw operation. The 
specific characteristics are as follows: 
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r^-M 


System Frequencies: 


Rotor Characteristics: 


Rotor Radius 16.67 ft 
Blade Chord (constant) 1.5 ft 
Coning Angle 3.5° 

Blade Twist 0.0° 

System Frequencies: 

1st Bending (fore-aft) 15.1 rad/s 
2nd Bending (fore-aft) 53.1 rad/s 
1st Side to Side 15.9 rad/s 


1st Bending (fore-aft) 2.75 rad/s 
2nd Bending (fore-aft) 12.8 rad/s 
1st Side to Side 2.9 rad/s 
2nd Side to Side 9,5 rad/s 

Aerodynamic Properties: 

Lift Curve Slope 5.73 
Drag Coefficient, Cp« .008 
Stall not Modeled 


Operating Conditions: 


Aerodynamic Properties: 

Lift Curve Slope 5.7 

Drag Coefficient, Cq^ .02 

Stall not Modeled 

Operating Conditions Used: 


Hind Velocity (mean) 

16. 

.63 

MPH 

Rotor Speed (mean) 

73" 

.35 

RPM 

Pitch Setting (to ZLL) 

3. 

.0° 


Turbulence Length Scale 

300 


ft 

Rms Turbulent Intensity 

2. 

,03 

ft/s 

Approximate Output 

6 


kW 


Selected analysis results for this situation are 
shown in the power spectral density plots of 
Figures 4 through 6. The figures clearly illus- 
trate the difference in response for the two wake 
models. For forces and moments which are highly 
dependent on the streamwise velocity such as 
thrust, the equilibrium wake assumptions give a 
larger response, as shown in Figure 4. This 
figure also indicates a significant response at the 
two system fore-aft bending frequencies. Figure 
5 shov/s the yaw response, while Figure 6 gives 
the pitching moment response. These two figures 
indicate the degree of coupling between pitch and 
yaw for this free -yaw turbine. Notice that as 
the yaw response increases, there is a correspond- 
ing increase in the pitching moment. Figure 6 
also shows a small response peak at 3Q, which is 
the result of the sinSfit and cos3Qt In the and 
Yy' inputs. 

The second wind system to be analyzed in this 
study was a large turbine called the Mod-G. The 
Mod-G was 2.5 MW turbine with a three bladed rotor 
located upwind of the tower, and was designed for 
fixed-yaw operation. The specific characteristics 
of this system are: 

Mod-G 

Rotor Characteristics: 

Rotor Radius 

Blade Chord (linear taper) 


Coning Angle 
Blade Twist (linear) 


150 ft 
7.74 ft 
at hub to 
3.15 ft 
at tip 
4° 

8 ° 


Wind Velocity 

20 


MPH 

Rotor Speed 

17. 

,5 

RPM 

Pitch Setting at tip 

-6. 

,2° 


Turbulence Length Scale 

500 


ft 

Rms turbulent intensity 

2. 

.44 

ft/; 

Approximate Power Output 

1, 

,1 

MW 


Figures 7, 8 and 9 shows $ome typical results for 
the Mod-G. As was the case in the previous plots, 
the system frequencies are easily identified. 

The primary objective of this work was to identify 
the features of turbulence which are most important 
In wind turbine design. In an effort to focus on 
these key features, the response at specific sys- 
tem frequencies was broken down into fractional 
contributions from each turbulence input. The 
most significant results of these calculations 
are tabulated in Tables 1 and 2. 

From these results it would appear that the most 
Important inputs are the three longitudinal turbu- 
lence terms, and in some instances, the two 1n- 
plane shear terms which have a oMfecti ye fre- 
quency of 2Q. An alternate means of presenting 
this same Information Is to plot power spectral 
density curves for the various outputs using only 
the turbulence input Vy, and then adding the two 
gradients Vy x and Vy,z. This has been done in 
Figures 10 through 15 for the outputs previously 
presented. As can be seen from these figures, 
a major contribution to the machine excitation 
is lost if the turbulence gradients are not in- 
cluded. However, neglecting the yy, and Cy, appears 
to have only a local influence around the frequen- 
cy 3Q. 

Conclusions and Recommendations 


On the basis of the work done in this study, the 
longitudinal turbulence input, Vy, and the two 
gradients, Vy,x and Vy,z are of equal importance 
when computing the dynamic response of wind sys- 
tems, and these three Inputs together comprise 
the major excitation source for horizontal axis 
wind turbines. Because of the simplifying assump' 
tions and approximations used in this analysis, 
it is Imperative that the results and the tech- 
nique be validated with experimental data, prior 
to use for design. 
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PSD OF VAW ANGLil Of THRUST 



Frequency /V rad/s Frequency rad/s 


ifeD'fl POWER SPECTRAL DENSITY PLOT OF THRUST. F^, FIGURE 6. MoD'M POWER SPECTRAL DENSITY OF PITCH NOHENT. 


ERRATUM: The power spectral densities 
shown in Figures 4-15 are Incorrectly 
plotted. The ordinate, Psd, should be 
increased by exactly a factor of 10 to 
be correct. 



o 

a 

a_ 



Frequency ~ rad/s 


Figure 5, Mod-11 power spectral density of yaw angle. (5. 


Figure 7. [Iod-G power spectral density of thrust. F^. 
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PsD OF SIDE FORCE 




Figure 8. Mod*G power spectral density of side force, Figure 10. The effect of the gradients i 

THRUST FOR MOD-M USING THE EQUILIBRIUM WAKE. 


ERRATUM: The power spectral densities 
shown in Figures 4-15 are incorrectly 
plotted. The ordinate, Psd, should be 
increased by exactly a factor of 10 to 
be correct. 




Figure 9. Mod-G power spectral density of yaw moment, Figure 11. The effect of the gradients x ^ i 

ANGLE FOR flOD-M USING THE EQUILIBRIUM WAKE. 
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•y-nt'i-mj ^ |^XTTnTfT[“Tnf 



Frequency ^ rad/s 


The effect of the gradients ^ 2, ^ oh pitch 

:W£NT FOR i'IcD-H USING THE EQUILIBRIUM WAKE. 



Frequency ^ rad/s 

Figure 1^. The effect of the turbulence gradients 

ON the thrust for the Mod-6 using the equilibrium 

WAKE. 


ERRATUM: The power spectral densities 
shown in Figures 4-15 are incorrectly 
plotted. The ordinate, Psd, should be 
increased by exactly a factor of 10 to 
be correct. 




:gure 13. The effect of the gradients on side Figure 15. The effect of the gradients on yaw 

FORCE FOR the MoD- 6 USING THE EQUI LIBRIUM WAKE. MOMENTS FOR THE MoD-6 USING THE EQUILIBRIUM WAKE 



Table 1. Fractional response contributions of the turbulence Inputs for the Mod-M using the 
equilibrium wake. 


Response/ Input 

% 






other 

Frequency 0 







Side Force, Fx 

.0 

.96 

.0 

.0 

.0 

.04 

Thrust, Fy 

1.0 

.0 

.0 

.0 

.0 

.0 

Pitch Moment, M 

.0 

.82 

.12 

.0 

.0 

.06 


Frequency = 15.1 (1st Bending) 


Side Force, F^^ 

.33 

.07 

.56 

.0 

.0 

.04 

Thrust, Fy 

.75 

.06 

.18 

.0 

.0 

.01 

Pitch Moment, M 

.70 

.07 

.22 

.0 

.0 

.01 

Frequency = = 23 

Side Force, 

.0 

.02 

.69 

.14 

.13 

,02 

Thrust, Fy 

.32 

.08 

,52 

,04 

.04 

.0 

Pitch Moment, M 

.0 

.11 

.77 

.06 

.06 

.0 


Table 2. Fractional response contributions of 
equlibrium wake. 

the turbulence 

inputs for the Mod-G 

using the 

Response/ Input 


V.X 

.Vz 



Other 

Frequency - 0 

Side Force, Fv 

.0 

.06 

.92 

.0 

.0 

.02 

Thrust, Fy 

1.0 

.0 

,0 

.0 

.0 

.0 

Yaw Moment, 

.0 

.97 

,0 

.0 

.0 

.03 

Pitch Moment, 

.0 

.0 

.97 

.0 

.0 

.0 

Frequency = 2.76 (1st 
Side Force, Fx 

Bending) 

.90 

.02 

.07 

.0 

.0 

.01 

Thrust, Fy 

.77 

.0 

.22 

.0 

.0 

.01 

Yaw Moment, M^ 

.77 

.01 

.21 

.0 

.0 

.01 

Pitch Moment, Mx 

.76 

.0 

.23 

.0 

.0 

.01 

Frequency = 3n = S. 5 

Side Force, Fx 

.01 

.36 

.07 

.27 

.27 

.02 

Thrust, Fy 

.06 

.05 

.42 

.22 

.24 

.01 

Yaw Moment, M^ 

.01 

.29 

.0 

.35 

.34 

.01 

Pitch Moment, M 

.0 

.05 

.45 

.23 

.25 

.02 


96 



APPENDIX 


Governing Equations: 







/"’I 



1^ 








M 

m-ji 

0 

”l3 

0 ( 

} 

Li 



^11 

0 

^13 

S 4 

0 

u 


'll 

0 

M 

22 ® 

”24 * 

3 

V 



0 

S 2 

0 

0 

Sb 

V 


0 

"31 

0 

”33 

0 { 

3 




Si 

0 

^33 

S 4 

0 


► + 

0 

0 

M 

42 0 

M 44 < 

3 

X 



Si 

0 

S 3 

S4 

0 

X 


0 

0 

0 

0 

° ”5bJ 




0 

S 2 

0 

0 

Sb 

V J 


0 



n 


p,. 

0 

1 

'13 

^14 

^15 

0 

>^17 

fl8 “1 



M 


0 


F 22 

0 

0 

0 

S6 

0 

0 


S 9 

= 


Y 

+ 

^31 


0 


'33 

S4 
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Ss 
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^52 
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Inertia Matrix 

= '"ll ^ “"n + ^ ’ ^13 = ^31 = -^'"n 

'^22 = "■*22 * ""r ^ ’ ”24 = "’24 = ^42 ’’ ^33 = "'33 ^ 'zz 

H.* = m.. + I ; Mrc = I ; m = mass of rotor ; m = mass of nacelle 

44 44 XX 55 r r n 


I 


S 4 

?2 

0 

C\J 

> 

S 3 

^34 

'42 

S 3 

S 4 

) 

0 

0 





(e^^cos3at + y^^sinSnt) 




I'- J 


q = distance from C|_ tower to nacelle-rotor C.G. ; Ixx ^nd lyy = mass moment of inertia of nacelle- 
rotor system about x and y axes ; I^ = rotor effective spinning inertia ; m-jj = tower inertia 
coefficients of Eq. (4), where for a uniform cantilever tower. 


"99 m^/420 , ^22 =156 m^/420 , m 2 ^ = 22 m^L/420 
2 

^33 ’ ^44 ” ’ ^t ” ^ tower polar inertia 


uamping Matrix 

C,1 = 3f(B^+e^Fg)/2RQ ; = -3f (5B^+B^Fp/2fi ; = -3f(C,+5e^E^j)/2n 

C 22 = 3fFjj./RJJ ; Cgg = -3fE^/n ; = -3f(e^{F^+B*}+a{B^+BQF^})/2n 

C_3 =• 3fR{{F2+6QaB*}+a{aB^+B^Fp)/2fl ; = 3fR({E^a+6^C*}+a{C^+a6^jEQ})/2SHI^fl 

C41 = 3f({E^.B2cp.a-6^{C^+E^})/2Q ; C 43 = -C 34 ; = C 33 

C ^2 ^ 3fC^/0 ; Cgg = SfRB^/^^+Cg ; = Generator torque coefficient 


Stiffness Matrix 

" *^11 ’ ^13 ^ "3fG^/2 ; K ^4 = -3f3QH^/2 ; K 22 = l <22 ; ><24 = ^24 
K 33 = IC 33 . 3fR(3,G*+iGj/2 ; K 34 = 3fR{H^+3,a-H^)/2 ; K ,2 = K 24 
K43 - “^34 s 1^44 = ^44 **■ 3 fR{BQ 6 |'+aG ^)/2 

where k^j = tower structural stiffnesses from Eq. (3), and for a uniform cantilever tower kn = 3EI/L^ ; 
k22 = 12EI/L3 ; k24 = SEl/l^ ; k33 = GJ/L ; k44 = 4EI/L 
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Wind Input Matrix 

>^11 = 3f{eXFo)/2R« ’ ^13 = ’ ^14 = -3f6,F,/2fi 

F^5 = 3fG^/2£J ; = -3f6^j(C^ -E^ )/2f2 ; F^g = -3f (-B,+B qF^ )/2S2 

T = 3fD^ ; Fgg = 3fFp/Rfi ; Fgg = 3fE^/n ; Fgg = -3fF,B^/fi 
F 31 = -3f(6jF^+B*Hi{B^+62p^j)/2„ . = 3f ({E^+3^C|}+i6^{C^+E^})/2fi 

F34 = 3fR(F2+6j,5F^)/2fi ; F35 = -3fR(6QC*+iC^ )/2n 

F 37 = -3fR({E2-Bjc*}-iBjj{C^-E^})/2fi ; F 3 g = 3fR(SgCF2-B*}+i{-B,+BQF^ })/Zfi 

^41 " ''33 ’ ''43 " *''31 ■’ ^"44 " ''35 ’ ''45 " ■''34 ’ ''47 " ■''38 ’ '^48 " ''37 
Q = 3 fRA^ ; F52 = 3 fC^/n ; Fgg = - 3 fR Bg/n ; Fgg = - 3 RBgC 2 /n 

where f = i pa’Rc.(Rn)^. The single subscript capitalized coefficients through are integral aero- 
dynamic coefficients of the form 
R 

A = / A'(x) dx where n = 0 , 1,2 
n 

with A' through F' defined as given in Eq. ( 13 ) for the "frozen wake" or Eq. ( 18 ) for the "equilibrium 
wake". In addition, 6 '(x) = AB'U) and H'(x) = AE'(x), while the coefficients with stars are B^=Bn-bBn-l , 
C^=Cn-fiCn-l and GJ^=Gn-FiGn-l , and R=h/R, a=a/R. 
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QUESTIONS AND ANSWERS 


R.W. Thresher 


From: B.J. Young 

Q: Your results show a substantial excitation which is continuous down through zero 

frequency, while J.P. Sullivan’s results were zero at zero frequency. Any coiranent 
on differences? 

A: Professor Sullivan used the Davenport model for the horizontal component of wind 

turbulence which vanishes at zero frequency . \^e used the model suggested by Von 

Karman which is finite at zero frequency ; therefore , this excitation difference 
at low frequency is due to the turbulence models. 

From: K.H. Hohenemser 

Q: How should you expect the results to change for hinged (teetering) blades? 

A: I have not done the analysis ^ so I do not know exactly, I would expect the magni- 

tude of the forces and moments to decrease significantly , 

From: Anonymous 

Q: What is the effect of damping on the first and second blade bending power spectra 

spikes? 

A: The blade is rigid so there are no blade resonances. The bending resonances illus- 

trated in the plots are for the tower. The only damping in the model is aerodynamic 
damping^ hut if structural damping were added^ the response near resonance points 
should be reduced. 

From: R. E, Wilson 

Q: Do you plan to treat teetering rotors? 

A: I would like to add teetering to the models hut at this time, our sponsor has not 

indicated a strong interest in adding this additional degree of freedom. Perhaps 
after the model is validated this will be possible , 
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